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Efforts to find the magnetic monopole in real space have been made in cosmic 
rays and in accelerators, but up to now there is no firm evidence for its existence 
due to the very heavy mass ~ 10 16 GeV. However, we show that the magnetic 
monopole can appear in the crystal-momentum space of solids in the accessible 
low energy region (~ 0.1 — leV) in the context of the anomalous Hall effect. 
We report experimental results together with first-principles calculations on the 
ferromagnetic crystal SrRuC>3 that provide evidence for the magnetic monopole 
in the crystal-momentum space. 

In 1931, Dirac Ij postulated the existence of magnetic monopole (MM) searching for 
the symmetry between the electric and magnetic field in the laws of electromagnetism. The 
singularity of the vectorpotential is needed for this Dirac magnetic monopole. Theoretically, 
the MM was found (3, y| as the soliton solution to the equation of the non-Abelian gauge 
theory for the grand unification. However its energy is estimated to be extremely large 
~ 10 16 GeV, which makes its experimental observation difficult. In contrast to this MM in 
real space, one can consider its dual space, namely the crystal momentum (k-) space of 
solids, and the Berry phase connection |4( of Bloch wavef unctions. This magnetic monopole 
in momentum space is closely related to the physical phenomenon, i.e., the anomalous Hall 
effect (AHE) in ferromagnetic metals. 

The AHE is a phenomenon where the transverse resistivity p xy in ferromagnets contains 
the contribution due to the magnetization M in addition to the usual Hall effect. The 
conventional expression for p xy is 

p xy = R B + AnR s M (1) 

where B is the magnetic field, Rq is the usual Hall coefficient, and R s is the anomalous 
Hall coefficient. This expression implicitly assumes that the additional contribution is pro- 
portional to the magnetization, and is used as an experimental tool to measure the mag- 
netization as a function of temperature. Especially it is extensively used in the studies of 
ferromagnetic semiconductors with dilute magnetic impurities, which are the most promis- 
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le spintronics p|. However, the mechanism of AHE has been controversial 
The key issue is whether the effect is intrinsic or extrinsic and how 
to treat the impurity and/or phonon scatterings. Karplus and Luttinger jfl] was the first to 
propose the intrinsic mechanism of AHE, where the matrix elements of the current operators 
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are essential. Other theories ?| Q] attribute the AHE to the impurity scattering modified 
by the spin-orbit interaction, namely the skew scattering 7] and/or the side-jump mecha- 
nism jlO(. These extrinsic mechanisms are rather complicated and depends on the details 
of the impurities as well as the band structure. Nevertheless, all those conventional theories 
on □ for AHE derives Eq.l, as they are based on the perturbative expansion in the 
spin-orbit coupling (SOC) A and the magnetization M, i.e., R s oc A. 
Recently the geometrical 



nized by several authors 



11 



meanin g o f the AHE of the intrinsic origin [6] has been recog 
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141 ]. The transverse conductivity be written as 



the integral of the Berry phase curvature (gauge field) over the occupied electronic states 
in crystal momentum space (Eq.5 below). Magnetic monopole corresponds to the source or 
sink of the gauge field/curvature defined by this Berry phase connection. Therefore the AHE 
can be the direct fingerprint of the MM in crystal momentum space. Note here that the 
presence of the time-reversal symmetry results in a xy = p xy = in the d.c. limit from the 
very generic argument, and the group theoretical condition for the nonzero a xy is equivalent 
to that of finite ferromagnetic moment (see Supporting Online Material (SOM)). Therefore 
we need the ferromagnets to study a xy even though the Berry phase connection is more 
universal and exists even in nonmagnetic materials. 

In this paper, we show by detailed first principles band calculation combined with the 
transport, optical, and magnetic measurements, that the observed unconventional behavior 
of the anomalous Hall effect and Kerr rotation in metallic ferromagnet SrRuOs is of intrinsic 
origin and is determined by the MM in k-space. The conventional expression Eq.l is invali- 
dated by the experimental data presented below, showing the non-monotonous temperature 
dependence including even a sign change. 

SrRuOs with the perovskite structure is an itinerant (metallic) ferromagnet. There are 
4 t2 g electrons with the low spin configurations. The Ad orbitals of SrRuOs are relatively 
extended and the bandwidth is large compared with the Coulombic interaction. The rela- 
tivistic SOC is also large in Ad electrons because of the heavy mass (of the order of 0.3eV 
for Ru atom). High-quality single crystal is available with the residual resistivity of the 
order of 10/zf2cm. These aspects make this system an ideal candidate to observe the AHE 
due to the k-space gauge field. Stoichiometric SrRu0 3 (bulk and thin film) and Ca-doped 
(Sr .8Cao.2Ru0 3 thin film) single crystals were prepared (see SOM for the details of sample 
preparation). The magnetization curve (Fig.lA) of SrRuOs film is quite similar to that of 
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bulk single crystal, except the Curie Temperature T c (~150K) is slightly lower than that of 
bulk (~160K), due to the strain effects. On the other hand, the isovalent Ca-doping sup- 
presses the T c and magnetization M dramatically. All of the samples were used for the trans- 
port measurement in the d.c. limit. As seen from Fig.lC, the p xy changes non-monotonously 
with temperature including even a sign change. Such behavior is far beyond the expectation 
based on the conventional expression Eq.l. In addition to the transport measurement, the 
frequency (u) dependent conductivities (Fig. 2) were measured for SrRuOs film by optical 
method (see SOM). Except the strong structures around 3.0eV, which are mostly due to 
the charge transfer from 0-2p to Ru-4<i, sharp structures are also observed for both the real 
and the imaginary part of a xy (u) below 0.5eV. Those low energy sharp structures cannot 
be fitted by extended Drude analysis. The lower the energy is, the stronger the deviations 
from fitting are. We will show, by combination with the first-principles calculations, that 
those unconventional behavior actually originate from the singular behavior of MM in the 
momentum space. 

Now we turn to some details of the theoretical analysis. Berry phase is the quantal 
jhase acquired by the wavefunction associated with the adiabatic change of the Hamiltonian 
4,^|. Let \n(a) > be the ra-th eigenstate of the Hamiltonian H(a) with a = (ai, ....,a m ) 
being the set of parameters, the Berry's connection is the overlap of the two wavefucntions 
infinitesimally separated in a-space, i.e., 



< n(a)\n(a + Aa) >= 1 + Aa- < n(a)\V a \n(a) >= exp 
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where the vector potential a n (a) is defined by a n (a) = i < n(a)\V a \n(a) >. Although 
the concept of Berry phase has broad applications in physics |15| . its relevance to the band 
structure in solids has been recognized only recently and in limited situations such as the 
quantum Hall effect under strong magnetic field !16 and the calculation of electronic polar- 

nn 

ization in ferroelectrics [13, IXSI] ■ In this case the parameter a is the crystal momentum k. 
For the Bloch wavefunction ^ n ^{r) = ^ tkr u n ^(f) with n denoting the band index and u n ^ 
being the periodic part, the vector potential for the Berry phase a nfi (k) is 

d 

a>n„(k)=i<u n% \—u n] ;>. (3) 

With this vector potential, the gauge covariant position operator x M for the wavepacket 
made out of the band n is given by = id^ — a nA1 (k). Therefore the commutation relation 



of x M 's includes the gauge field F^ v = d kli a nv - d ku a nil as 

%J\ J-f/ii/' (4) 

which leads to the additional (anomalous) velocity —i[Xfj,,V(x)] = —F^dV(x)/dx u being 
transverse to the electric field E u = —dV{x)/dx u . Therefore the transverse conductivity o xy 
is given by sum of this anomalous velocity over the occupied states as |3] 

cr xy = J2 n F( £ n( k -))b z (k). (5) 

n,k 

where b z (k) = F xy (k) and Uf{s) = l/^e^ 6- ^ + 1) ( (3: inverse temperature, fi: chemical 
potential ) is the Fermi distribution function. Hence the behavior of gauge field b z (k) in 
fc-space [19] determines that of a xy . One might imagine that it is a slowly varying function 
of k, but it is not the case. Fig. 3B is the calculated result for b z (k) in the real system 
SrRuOs. It has a very sharp peak near T-point and also ridges along the diagonals. The 
origin for this sharp structure is the (near) degeneracy and/or the band crossing, which act 
as MM. Consider the general case where two band Hamiltonian matrix H(k) at k can be 
written as H(k) = X^=o, 1,2,3 //^(k)cr M where 0X2,3 are the Pauli matrices and a is the unit 
matrix. We can consider the mapping from k to the vector f (k) = (fi(k), /2(k), /^(k)), and 
the contribution to a xy from the neighborhood of this degeneracy region is given by the solid 
angle dQf for the infinitesimal dk x dk y integrated over k. Therefore the gauge flux near the 
MM, namely the degeneracy point f = is singularly enhanced (as shown in Fig. 3). (See 
SOM for more details). 

We studied the behavior of o xy by the first-principles calculations (see SOM for meth- 
ods). The calculated density of states (DOS) is not so different between the cases with and 
without SOC (Fig.4A), while the a xy should be very sensitive to the Bloch wavefunctions 
and depends on the Fermi level position and the spin-splitting (magnetization) significantly, 
as predicted by above discussions. The behavior of function of the Fermi-level 

position was obtained by using a small broadening parameter for the lifetime 5 (=70meV) 
(Fig.4B). By shifting the Fermi level, not only the absolute value but also the sign of a xy 
is found to change. The sharp and spiky structures are just the natural results of singular 
behavior of MM (Fig.3). For the case without any shift of Fermi level, we obtain the value 
of a xy = —60 f2 _1 cm _1 , which has the same sign as and is comparable with the experimental 
value (~ — lOOf^cm -1 ). Such a spiky behavior should be also reflected in the ui dependent 
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a xy , especially for the low energy range with longer lifetime, while be suppressed at higher 
activation energy with shorter lifetime. As shown in Fig.2, for the u dependence of optical 
conductivity, the high energy part (>0.5eV), which is dominated by the p — d charge transi- 
tion peak, is usual, and can be well reproduced by our calculations. While the observed peak 
structure of a xy (u) below 0.5eV is the clear demonstration of the predicted spiky behavior. 
The spectra below 0.2eV is not measured due to the technical difficulty, but it is clear that 
even sharper structure should be there because the d.c. limit Ke(a xy ) ~ — lOOf^cm -1 has 
the opposite sign (the lm(a xy ) at d.c. limit should go back to zero). Such a low energy 
behavior is well represented by our calculations, providing further evidence for the existence 
of magnetic monopoles. 

Now it is straightforward to understand the results of transport measurement for a xy . 
Here we attribute the temperature (T) dependence of a xy to that of the magnetization 
M(T). As the results of /c-space integration over occupied states, the calculated a xy is non- 
monotonous as a function of magnetization (Fig. ID). With the reduction of spin-splitting, 
the calculated a xy , after the initial increase, decreases sharply, then increases and changes 
sign becoming positive, and finally goes down again, capturing the basic features of the ex- 
perimental results. Even more surprisingly, by converting the measured p xy versus T curves 
shown in Fig.lC into the a xy versus M curves shown in Fig. ID, now they all follow the 
same trend, and match with our calculations as long as the experimental data is available. 
Those curves are measured for different samples (with different saturation moments), they 
all follow the same rule qualitatively, and could be simply explained by the reduction of mag- 
netization |2]J (see Fig.lA). Note, however, that the comparison between the experiments 
and calculations should be semi-quantitative, because the results are sensitive to the lattice 
structures. The calculated a xy for the fictitious Cubic structure shows a strong deviation 
from that obtained for orthorhombic structure, and it changes the sign to be positive at low 
temperature (large M). Therefore more accurate information on the structure is needed to 
obtain the quantitative result. However, such a sensitivity does not affect our main results, 
i.e, the non-monotonous behavior of a xy . Even the calculations for Cubic structure show 
such a behavior, and may be used as a guide of possible deviation. 

The results and analysis presented here should stimulate and urge the reconsideration 
of the electronic states in magnetic materials from the very fundamental viewpoint. For 
example, the magnetic monopole is accompanied by the singularity of the vector potential, 
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ring 



i.e., Dirac string As shown by Wu-Yang 



22j this means that more than two overlapping 



regions have to be introduced, in each of which the gauge of the wavefunction is defined 
smoothly. This means that one cannot define the phase of the Bloch wavefunctions in a 
single gauge choice when the magnetic monopole is present in the crystal momentum space. 
This leads to some nontrivial consequences such as the vortex in the superconducting order 
parameter as a function of k j^J , and many others are left for future studies. 
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FIG. 1: Measured temperature dependence of the (A) magnetization M, (B) longitudinal resistivity 
p xx , and (C) transverse resistivity p xy for the single crystal and thin film of SrRuOs, as well as for 
the Ca doped Sro.sCao^RuOs thin film. The corresponding transverse conductivity a xy is shown in 
(D) as a function of the magnetization, together with the results of first principle calculations for 



cubic and orthorhombic structures 



201 ] . In our calculations, the change of magnetization is taken 



into account by the rigid splitting of up and down spin bands. As the transverse conductivity 
should vanish with M at high temperatures, the calculated a xy is multiplied by the additional 
M/Mq (Mq = 1.5/ie) factor, which does not affect its behavior except at the very vicinity of Tc. 
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FIG. 2: Calculated (left panels) and measured (right panels) longitudinal (a xx ) and transverse {<J xy ) 
optical conductivity of SrRu03 film. The measurements were performed at low temperature (10K). 
The calculations were done for both the orthorhombic single crystal structure and the hypothetic 
cubic structure by keeping the average Ru-0 bond length. The experimental a xy are shown by 
multiplying a factor of 4. The quantitative comparison between the experiments and calculations 
about the absolute value of a xy would require more accurate structure information (see text part). 
Nevertheless, the clear peak structures for the low energy a xy is the demonstration of monopoles 
associated with b z (k). 10 
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FIG. 3: (A) Geometrical meaning of the contribution a xy when the two bands are nearly de- 
generate. The two-dimensional Hamiltonian matrix H{k) can be generally written as H(k) = 
S^=o 123 fni^) a fi where (71,2,3 are the Pauli matrices and o"o is the unit matrix. By mapping from 
k to the vector f(k) = (f i(k) , f 2(h) , f 3(h)) = f(k)(cosippsm9p,shnppsva9p,cos9p), the contribu- 
tion to a xy from the neighborhood of degeneracy region can be given by the /-space solid angle 
(Klj- = [d(9 ptp j)/ d(k x ,k y )] sin 6 jdk x dk y = dip ^-sm 9 pl9 ^ ior the infinitesimal dk x dk y integrated 
over k. The solid angle corresponds to the flux from the monopole at / = 0. (See SOM). (B) 
Calculated flux distribution in /c-space for t% g bands as a function of (k x ,k y ) with k z being fixed 
at for SrRuOs with Cubic structure. The sharp peak around k x = k y = and the ridges along 
k x = ±k y are due to the near degeneracy of d yz and d zx bands because of the symmetry reasons 
as explained in the online supporting text. 
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FIG. 4: The calculated (A) density of states (DOS) and (B) a xy as functions of Fermi level position 
for the orthorhombic structure of single crystal SrRuOs. The Fermi level is shifted rigidly relative 
to the converged solution, which is specified as zero point here. The sharp and spiky structure of 
a xy is the demonstration of the singular behavior of magnetic monopoles. 
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SUPPORTING ONLINE MATERIAL 
Materials and Experimental Methods 

Both stoichiometric SrRu0 3 and Ca-doped (Sr . 8 Ca .2RuO 3 ) thin (500A) single crystal 
epitaxial films were grown on the (001) surface of SrTi0 3 (STO) single crystal substrate by 
pulsed laser deposition (PLD) employing KrF laser pulses (lOOmJ) focused on the polycrys- 
talline target, while SrRu0 3 bulk single crystal were prepared using a flux method (SrCl2 
flux). The quality of those samples were confirmed by the X-ray diffraction. The bulk sin- 
gle crystal is orthorhombic, with a=3.91lA, 6=3.936A, c=3.922A. The films are coherently 
strained by the SrTi0 3 substrate, yielding a tetragonal distortion in the [001] direction. As 
the results, the out-of-plane lattice constants of films are elongated (c=3.950A for SrRu0 3 
film and c=3.93lA for Ca-doped film), and the magnetic easy axis is perpendicular to the 
film plane. 

For the transport measurements of films, they were patterned in Hall bar geometry us- 
ing conventional photo-lithography and Ar ion dry-etching. The Hall resistivity pu was 
measured together with the longitudinal resistivity p xx as a function of temperature under 
applied magnetic field. The anomalous resistivity p xy was determined after subtraction of 
the ordinary Hall contribution from the measured pn, and the transverse conductivity a xy 
was estimated as —p xy / p 2 xx - In addition to the transport measurement, the frequency (o>) de- 
pendent conductivities were measured for SrRu0 3 film by optical method. The longitudinal 
part of the optical conductivity a xx (uj) was obtained from Kramers-Kronig transformation 
of normal reflectivity from 0.08eV to 40eV, while the transverse part cr xy ((v) was deduced 
from magneto-optical Kerr spectra. The spectra were measured by a polarization modu- 
lation method. For the energy region of 0.7 - 4 eV, ordinary measurement system with 
grating monochromator and photo-elastic modulator (CaF 2 window, modulation frequency 
of 56 kHz) was used for the measurement. For the energy region below 0.8 eV, Fourier 
transform infrared (FT-IR) spectrometer with rapid scan mode (scan speed 0.16 cm/s) was 
utilized for the measurement. The light from the FT-IR interferometer was, at first, polar- 
ized by wire-grid polarizer (BaF 2 substrate). Then the polarization was modulated by ZnSe 
photo-elastic modulator (modulation frequency of 50 kHz) and focused by BaF 2 lens (focal 
length of 150 mm). The light reflected from the sample was detected by photovoltaic type 
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HgCdTe detector. By this system we can measure magneto-optical spectra down to 0.2 eV 
with accuracy higher than 0.01 degree. 



Some Details of First-Principles Calculations 

First-principles calculations of a xy are quite challenging, requiring the combination of sev- 
eral modern calculation techniques and powerful supercomputer sources. The plane-wave 
pseudopotential calculations jl| were performed based on the local spin density approxima- 
tion (LSD A), that give good descriptions for the electronic and magnetic properties of this 
compound P, iij]. The spin-orbital coupling (SOC) is treated self-consistently by using the 
relativistic fully separable pseudopotentials |4j in the framework of non-collinear magnetism 
formalism. The inter-band optical transitions are calculated from the converged Kohn-Sham 
eigea states by using the Kubo fonna.a QQ The analytical tetrahedron raethe, Q has 
been used for the accurate fc-space integration, and the convergence has been checked care- 
fully. The finite life-time broadening 5 has been estimated from the experimental residual 
resistivity and extended Drude analysis of a xx . The validity of our calculation techniques 
has been well demonstrated in the explanation of anisotropic optical data in Ca 2 Ru0 4 ^]. 

a xy and Magnetic Monopole in Momentum Space 

Here we present more details on the contribution to the transverse conductivity a xy 
from the k region where two bands are nearly degenerate. Enhancement of the gauge field 
b z (k) occurs when more than two bands are close energetically, which corresponds to the 
magnetic monopole and shows its fingerprint in anomalous Hall effect (AHE) and Kerr 
rotation as described below jjj. Consider the general case where two band Hamiltonian 
matrix H(k) at k can be written as H(k) = Z)u=oi,2 3 fnitycru where cri,2,3 are the Pauli 
matrices and <To is the unit matrix. Then we can consider the mapping from k to the vector 
f(k) = (f i(k) , f 2(h) , f 3(h)) = f(k)(cos </?j*sin 9? : sin <^?sin 9 ? : cos 9^) as shown in Fig. 3A of 
the main text. Then H{k) can be easily diagonalized to obtain the two eigenvalues £±(k) = 
fo(k) ± f(k). Calculating Eq.5 of the main text in this case, we obtain the contribution to 
a xy from these two bands as 



e 



2 



8nh 



<y = 5-r d*k[n F {e-(k)) - n F {e + {k))] 
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, dk x dk v dk v dk x J * 



e 2 



f dk z dn f [n F (e4k))-n F (s + (k))]. (SI) 

Here dQj- = [d(9f, (pf)/d(k x , k y )] sa\dfdk x dk y = dtp rsin Qfddj is the /-space solid angle, 
which is the integral of the gauge field 

b(f) = ±-t- 
l/l 3 

due to the monopole at / = over the infinitesimal surface in /-space corresponding to 
the small square dk x dk y in fc-space (Fig.3A of the main text). Therefore o"^y bands again has 
the geometrical meaning in / space. This gauge field strongly depends on k in the (near) 
degeneracy case, i.e., when f(k) is near the monopole. 

There are two cases for the (near) degenerate bands. One is the accidental degeneracy 



ll| where the three equations fi(k) = f2{k) = fz{k) = are satisfied at k = k . Near 
this accidental band crossing, one can expand as f a (k) = J2b a ab(kb — k ob ) where a, b = 1, 2, 3. 
In this case the /-space can be identified with the fc-space, and the gauge field distribution 
in fc-space around k = ko is similar to Eq.Sl above replacing / by k — k$. Although there 
occurs no singularity in a xy at /x = £±(&o) due to the cancellation between the positive and 
negative k z — k z0 , the nonlinear dependence of f a (k) gives rise to the strong /i-dependence 
of a xy slightly away from this energy. The other class of (near) degeneracy is due to the 
symmetry, where the fc-group has the irreducible representation with the dimensions more 
than 2. As an example, one can consider the simplest tight binding model of t2 9 orbitals 
d yz and d zx on the cubic perovskite structure, which is relevant to the SrRuC>3 discussed 
in the main text. In this case the H(k) for these two bands with up-spin is given by 
fo(k) = -2ticosk z - t^cos k x + cos k y ) , j\(k) = 2t 2 smk x smk y , f 2 (k) = -AM, f 3 (k) = 
—ti(cosk x — cosky), where ti,t% are the effective intra- and inter-orbital transfer integrals 
respectively, A is the spin-orbit coupling (SOC) constant, and M is the magnetization. 
When AM = 0, there occurs the degeneracy along the line k = (0, 0,k z ). Furthermore 
when £2 = 0, there occurs the degeneracy along the plane k x = ±k y . Considering the case 
t\ » ti >> AM, which is relevant to SrRu03, the gauge field b z (k) with k z = has a 
largest peak at k x = k y = and is enhanced along the lines k x = ±k y , which is actually 
seen in the realistic calculation for SrRu03 shown in Fig.3B of the main text. In this case 
there occurs no cancellation of b nz (k) from the integral over k z . 

15 



Although the discussion above is applicable to any (near) degeneracy, the (singular) gauge 

— * 

field from different band crossings in /-space cancel out in the presence of the time-reversal 
symmetry and/or the absence of the SOC. The former prohibits the finite a xy . In the absence 
of SOC, up-spin and down-spin bands are decoupled, and each of them can be represented by 
the Hamiltonian matrix without the time-reversal symmetry breaking. In ferromagnets with 

the SOC, the singular behavior of a xy by changing some parameters such as the chemical 

— * — * 

potential and the magnetization is the fingerprint of the monopoles in k- and/or /-space. 
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